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ABSTRACT 


The authors recently reported that maximum atmospheric densities measured 
by the Explorer XIX and XXIV satellites occurred in the early afternoon in the 
winter hemisphere, that is, the hemisphere opposite the equator from the sun. 
They referred to this density peak as the winter bulge. Previously, it was gen- 
erally assumed that maximum densities at these altitudes (550 km to 750 km), 
where atomic oxygen and helium are the primary constituents, would occur in the 
early afternoon in the summer hemisphere. In the present paper, drag data from 
the Explorer IX, XIX, and XXIV satellites have been used to obtain a more com- 
plete picture of the winter bulge. These data cover the altitude range 590 km 
to 880 km and the time interval between 1961 and 1966. A correlation is found 
between the bulge latitude and the inferred helium to atomic oxygen ratio. This 
correlation indicates that the winter bulge is the result of a high latitude 
winter helium bulge. Superimposed upon this helium bulge is a near-equatorial 
summer oxygen bulge. The winter helium bulge occurs in the early afternoon and 
migrates latitudinally in the opposite sense of the sun from high latitudes in 
the hemisphere opposite the sun at the solstices to near the equator at the 
equinoxes . A study of drag data from the Echo IX satellite confirms the high 
latitude winter helium bulge. The abundance of helium at high latitudes in the 
winter is not attributed to higher exospheric temperatures but to a seasonal 
variation in helium concentration at 120 km which may be related to the level 
at which diffusive separation of helium begins. 


THE DISTRIBUTION OF HELIUM AND ATOMIC OXYGEN 


IN THE LOWER EXOSPHERE 
By Gerald M. Keating and Edwin J. Prior 


INTRODUCTION 


Peak atmospheric densities occur at satellite altitudes at approximately 
2:00 every afternoon. This region of peak densities is termed the diurnal bulge. 
Although the local solar time of the diurnal bulge is generally agreed upon, the 
latitudinal migrations of the bulge center are still in question. Originally, 
peak atmospheric densities and exospheric temperatures were assumed to follow 
the latitudinal migrations of the sun (ref. l). More recently, this picture has 
changed after analysis of drag data from the Air Density Explorer satellites. 
Jacchia and Slowey (ref. 2) found they could obtain a better prediction rule by 
assuming the diurnal bulge remained on the equator rather than following the 
latitude of the sun. Keating and Prior (ref. 3 ) analyzed data from the latest 
two Air Density Explorers, Explorer XIX and XXIV, during a period of low solar 
activity between the altitudes of 550 km and 750 km. They found the diurnal 
bulge occurred in the winter hemisphere, that is, the hemisphere opposite the 
equator from the sun. The authors referred to this density peak as a "winter 
bulge." Using the winter bulge gives better agreement with observations than 
the equatorial bulge of reference 2 or the summer bulge of reference 1. The 
present study was undertaken to determine over what altitude range and/or range 
of solar activity a winter bulge existed. Densities inferred from the atmo- 
spheric drag upon satellites Explorer IX, XIX, XXIV, and Echo II were used for 
determining the latitudinal variations of the diurnal bulge over the years 1961 
to 1966 and over the altitude ranges of 390 km to 1130 km. A shift of bulge 
position with altitude appears to be related to the diffusive separation of atmo- 
spheric constituents in the exosphere. As a result, the study provided new 
insight into the distribution of the major constituents of the lower exosphere, 
atomic oxygen and helium. The high-altitude winter bulge appears to be the 
result of an increased helium concentration in the winter hemisphere while the 
observed summer bulge of . atomic oxygen may coincide with high exospheric 
temperatures . 


DENSITY DATA 


Four low mass-area ratio satellites were used in this study of the latitu- 
dinal variations of the diurnal bulge in the exosphere. These satellites were 
chosen due to their sensitivity to the drag environment. Explorer IX, XIX, and 
XXIV were each 12-foot-diameter ( 3 .66 m) Air Density Explorer satellites which 
weighed less than 10 kilograms. Echo II is a 135-foot -diameter (4l.l m) sphere 
which weighed 256 kilograms. Explorer IX was launched in February 1961 and 
reentered in 196 k. Densities near its perigee point were deduced from its energy 
decay. Its perigee point precessed in such a fashion that it was in near reso- 
nance with the relative motion of the sun. As a result, densities were measured 
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near the hour angle of the diurnal bulge throughout the satellite lifetime. 

Thus, it was a valuable tool for measuring variations in bulge position. In this 
study, Explorer IX measurements between the latitudes of ±39° were made over the 
altitude interval 390 3m to 880 3m. The Explorer XIX and XXIV satellites were 
placed into near-polar orbits. As a result, although their perigees do not 
always remain close to the diurnal bulge, they can adequately measure latitudinal 
variations of atmospheric density in the exosphere. ' The Explorer XIX and XXIV 
data studied include recent data not included in the study of reference 3- 
Explorer XIX data cover the latitudes ±79° and the altitude interval 630 3m to 
750 3m. Explorer XXIV data cover the latitudes ±81° and the altitude interval 
550 3m to 660 3m. The Echo II density data are taken from reference 4. It has 
been normalized to an altitude of 1130 3m, and covers the latitude range ±82°. 

The Echo II satellite orbit is sometimes so circular that the position of the 
density measurement is difficult to pinpoint. Consequently, the period of mea- 
surements when the eccentricity was less than 0.015 has been deleted. This 
covers the interval July 28, 1964, through January 5> 1965* 

The methods developed for the determination of densities from the Air Den- 
sity Explorer satellites have been described in detail by Keating, Mullins , and 
others (refs. 5 and 6). The effects of atmospheric drag are deduced by sub- 
tracting energy changes due to radiation forces from the total energy decay. 
Radiation forces were evaluated from reflectance measurements of the satellite 
material (refs. 6 and 7)» The assumed drag coefficient increases with perigee 
altitude. The density at perigee is evaluated by matching the observed energy 
decay due to drag with the energy losses obtained when the satellite is moving 
through the rotating spring/fall atmosphere of reference 8. This model atmo- 
sphere has a different density profile for each exospheric temperature, and is 
nearly identical with the static diffusion models of reference 1. In order to 
minimize errors in atmospheric density due to inaccuracies in the assumed atmo- 
spheric model, the density is then evaluated at an altitude between l/2 and 

\[S/2. density scale heights above perigee, depending on the value of the density 
scale height and its rate of variation (ref. 6). This is similar to the tech- 
nique employed in reference 9* 


THE DIURNAL BULGE 


The summer bulge model of Jacchia (ref. l), the equatorial bulge of Jacchia 
and Slowey (ref. 2), and the winter bulge of Keating and Prior (ref. 3) can all 
be obtained by varying the parameter B in the following expression given by 
Keating and Prior in equation (4) of reference 10: 


T = T 0 


1 + 0.28 sin m 0 T + 0.28(cos m V - sin m 9 ’ )cos^*5 ( X 
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where 

T)' - Ij&p - B6 g | , deg 
9' = ^J&p + B8gJ , deg 

t = H - 45° + 12° sin(H + 4^°), deg -si < r < jt 

H = hour angle of sun, deg 

5p = declination or latitude of density point, deg 
Sg = declination or latitude of sun, deg 
T 0 = minimum exospheric temperature, °K 

5p = BSg = declination or latitude of diurnal bulge center, deg 
B = constant 
m = coefficient 

T = exospheric temperature at H,Sp before geomagnetic activity correction, °K 

After calculating the exospheric temperature, the atmospheric density can 
be determined at any altitude by entering the spring/fall tables in the U.S. 
Standard Atmosphere Supplements, 19 66 (ref. 8) or the nearly equivalent tables 
of reference 1. The models give density as a function of altitude for any exo- 
spheric temperature. Actually, exospheric temperatures are deduced from the 
models and are therefore no more valid than the assumptions made to generate the 
models such as invariant boundary conditions at 120 km. However, inferred exo- 
spheric temperature is a useful parameter for predicting densities over a wide 
range of altitudes, and gives some indication of the range of temperatures in 
the lower exosphere. 

A study of the parameter B of equation (l) gives insight into the nature 
of the latitudinal migrations of the* diurnal bulge. Since the latitude of the 
bulge center is given by B times the latitude of the sun, a negative B indi- 
cates a diurnal bulge remaining in the winter hemisphere while a positive B 
indicates a diurnal bulge in the summer hemisphere. Shown in figure 1 is the 
variation of bulge position as a function of latitude and time of year for dif- 
ferent values of B. For each value of B, the local solar time of the bulge is 
constrained in equation (l) to approximately 1400. The diurnal bulge of refer- 
ence 1, which follows the latitudinal migrations of the sun, corresponds to 
B = +1 and m = 2.5- The equatorial bulge of reference 2 corresponds to B = 0 
and m = 1.5* The winter bulge of reference which gives even lower residuals 
than the first two models, corresponds to B = -1 and&f# = 1. The constant m 
indicates how rapidly density drops off latitudinally from the bulge center. A 
low value of m such as m = 1 indicates atmospheric density in the daytime 
drops more slowly from the bulge center with latitude than with hour angle. 
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In previous studies , a constant value of B has been assumed over the 
entire range of satellite data. It was therefore undertaken to determine 
whether the assumption of constant B was valid; and if B did vary, what 
parameters would affect it. Density data from the Air Density Explorer satel- 
lites and Echo II satellites were divided in such a fashion as to give represen- 
tative samples of latitudinal variations. For example, Explorer IX data were 
divided into approximately 5-month intervals corresponding to two revolutions 
of the perigee point. Thus, the range of latitudes was sampled four times in 
each interval. 

Densities were predicted for each observation, taking into account the 
effects of geomagnetic activity, solar activity, and the semiannual variation 
given in the U.S. Standard Atmosphere Supplements, 1966 (ref. 8). As a result, 
the minimum exospheric temperature on any date was assumed to be given by the 
following expression: 


T q = 362° + [3.60° F 10 7 + 1.8°(f 


- F 

10. T 10.7; 


+ 1.0 L 


0.57 + 0.14 sin ) 


F i0 7 sin kit 


a - 59 ) 
i 365 / 
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where 

10.7-cm flux from sun divided by 10“ 22 watts/m 2 /cycles/sec 

10. 7- cm flux from sun averaged over three solar rotations divided by 
10" 22 watt s/m 2 / cycles/sec 

d = number of days counted from January 1 

The minimum exospheric temperature T 0 was then substituted into equation (l), 

m was set equal to 1, and B was set equal to values ranging from -3*9 to 3-9- 
The increase in temperature related to geomagnetic activity was determined using 
the following expression (ref. 8): 


F 10.7 " 
f 10.7 = 


T' = T + 1.0°(ap) + 100° jl - exp(-0.08 a p )] (3) 

where 

T' = corrected exospheric temperature, °K 
a p = geomagnetic planetary amplitude 

The spring/fall models of reference 8 were generated up to 1130 km in order 
to obtain predictions of Echo II densities. The models indicate that for the 
Echo II satellite, helium drag is generally 10 times as high as hydrogen drag. 
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For each observed, density of the four satellites, a corresponding predicted 
"density was obtained using equations (l), (2), and ( 3 ), and the spring/fall 
models of reference 8. For each value of B in equation (l), a standard devi- 
ation between predicted and observed densities was obtained for each interval of 
data points. The values of B giving the lowest standard deviation for each 
data set are tabulated in table 1. 

The first two columns of table 1 give the satellite and data interval of 
each data set studied. The second two columns give the average altitude and 
average inferred temperature of each set. The last column, l°gio[” n (He)/n(0)J , 

is determined knowing the average altitude and average temperature of each data 
set and using the models of reference 8. The variation in the helium-to-atomic- 
oxygen ratio is due to the diffusive separation of constitutents in the exo- 
sphere. Standard deviations of logqo of the atmospheric density for individual 

data sets averaged about 0.08. Mean densities were generally slightly greater 
than the models above 600 km, and less than the models below 600 km. 

Shown in figure 2 is the variation of B with log^Q [n(He)/n(o)J corre- 
sponding to the tabulation in table 1. On the left of the figure, atomic oxygen 
is the primary constituent, while on the right, helium is the primary constituent. 
Above logioQ^He )/n (0 )1 of 0.4, all 11 data points give negative B and there- 
fore indicate a winter density bulge. Each of these data points represents many 
months of satellite data. The remaining four data points, in an atmosphere pri- 
marily composed of atomic oxygen, indicate a summer density bulge . There is a 
general trend for the value of B to shift from positive in an atomic oxygen 
atmosphere to negative in a helium atmosphere. Remembering that positive B 
indicates a summer bulge and negative B indicates a winter bulge, it is appar- 
ent that the observed diurnal bulge shifts hemispheres depending upon which con- 
stituent is predominant. It is also apparent that a constant B model such as 
the equatorial bulge of reference 2 (B - 0) does not adequately describe the 
location of the diurnal bulge center in the lower exosphere. 


The shift in bulge latitude can be understood by envisioning two bulges 
existing simultaneously, a winter helium bulge occurring at high latitudes in 
the winter hemisphere, and a summer oxygen bulge occurring at low latitudes in 
the summer hemisphere. This concept is pictured schematically in figure 3 • 

This diagram is not intended to show the shape or the altitudinal variations in 
the bulges. The sun is shown near summer solstice in the Northern Hemisphere. 
Atomic oxygen is shown to peak at low latitudes on the sun side of the equator. 
Helium peaks at high latitudes in the hemisphere opposite the sun. As the sun 
migrates latitudinally down to the equator, so do both bulges. Then the sun, 
followed by the atomic oxygen bulge, moves into the Southern Hemisphere and the 
helium bulge shifts into the Northern Hemisphere. With increasing altitude, 
helium replaces atomic oxygen as the primary constituent, and the region of 
maximum total atmospheric densities shifts from the summer hemisphere to the 
winter hemisphere. It also follows that with increasing solar activity and exo- 
spheric temperatures, atomic oxygen predominates to greater altitudes, and the 
summer bulge will be observed at greater altitudes. 
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Empirical expressions have "been developed for the atomic oxygen bulge and 
the helium bulge. These expressions are in terms of a virtual exospheric tem- 
perature which may in turn be converted to atmospheric density using the static 
diffusion models of reference 8. The helium bulge is considered to be the result 
of a seasonal variation in the boundary conditions given for 120 km. Thus, the 
exospheric temperatures for helium are considered virtual since they are based 
on invariant boundary conditions. The distribution of helium is predicted by 
setting the B of equation (1) equal to -3-0, and m equal to +1. The atomic 
oxygen distribution is predicted by setting B equal to +0.25, and m equal 
to +1 . 

A sample calculation is made to show how the number densities of helium 
and atomic oxygen can be predicted at any position in the lower exosphere and 
how exospheric densities can be calculated. Determine the number densities of 
atomic oxygen and helium for July 25, 1965? for an altitude of 700 km, an hour 

angle of the sun of 60°, and a latitude of 5° • First, calculate T Q from 
equation (2) setting F 10> y = 70, F 10 ^ = 75, and d = 205. It is found that 

T Q = 589° K. Next, calculate the virtual exospheric temperature of atomic oxygen 
from equations (l) and (3) setting T 0 = 589 0 K, H = 60°, 6p = ^5°, Sg = 20°, 

B = +0.25, m = 1, and ap = 6. Finally, determine the number density of atomic 

oxygen entering the tables of reference 8 at T' = 782° K and altitude equal to 
700 km. The number density of atomic oxygen equals 7-72 X 10^ cm-3. To deter- 
mine the virtual temperature of helium setting T 0 = 589° K, H = 60°, 5p = ^5°, 
5g = 20°, B = -3-0, and m = 1, in equation (l), it is found that T = 684° K. 
Setting ap = 6 in equation (3), T* = 728° K. Entering the tables of refer- 
ence 8 at T' = 728° K and altitude equal to 700 km, the number density of 
helium equals 6. 70 x 10^ cm ”3. Neglecting the minor contributions of other con- 
stituents, this gives a total density of 6.50 x 10-18 g/cm3. 


DISCUSSION 


Rather than a seasonally dependent unknown heat source near the boundary of 
the exosphere, it is expected that the winter bulge is the result of a seasonal 
variation of helium concentration in the lower thermosphere. This seasonal vari- 
ation may be related to a shift in the effective altitude at which diffusive 
separation of helium begins. A slight increase in this altitude, possibly 
related to increased turbulence , may sharply reduce helium concentrations in the 
upper thermosphere and exosphere (ref. 11). There is also an indication that 
total densities are higher and temperatures lower at 120 km in winter than in 
the summer (ref. 8). For a heavy constituent such as molecular nitrogen, this 
seasonal variation can be dwarfed at greater altitudes by the effects of high 
summer temperatures. For a lighter constituent such as helium, which is not as 
strongly affected by temperature, such a variation could extend into the exo- 
sphere. The low latitude summer oxygen bulge is expected to correspond to high 
exospheric temperatures. The oxygen bulge given by the model is closer to the 
equator than the subsolar point. This may result from the ion drag effects upon 
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the neutral atmosphere suggested by Jaechia and Slowey in reference 2 or from 
seasonal variations at 120 km indicated by Champion (refs, 8 and 12). 


CONCLUSIONS 


A study was performed to determine the latitudinal migrations of the diur- 
nal density bulge of neutral constituents in the lower exosphere. This study 
was based on satellite drag data in the altitude range 390 km. to 1130 km and 
the time period 1961 to 1966. The following may be concluded: 

1. At altitudes where atomic oxygen predominates, the diurnal bulge 
remains at low latitudes in the summer hemisphere, that is, hemisphere on sun 
side of equator. 

2. The diurnal bulge remains in the winter hemisphere at greater altitudes 
where helium replaces atomic oxygen as the major constituent. The winter bulge 
migrates up to high latitudes near the time of the solstices and down to the 
equator near the time of the equinoxes, 

3. The observed phenomena appear to result from the simultaneous existence 
of a high latitude winter helium bulge and a low latitude summer atomic oxygen 
bulge. As a result, a simple two-bulge model has been developed to describe the 
lower exosphere, and can be used to predict atmospheric densities. 

4. With increasing altitude, as helium replaces atomic oxygen as the pri- 
mary constituent, the diurnal bulge shifts from the summer hemisphere to the 
winter hemisphere. With increasing solar activity and exospheric temperature, 
the peak altitude where atomic oxygen predominates will increase, and it is 
expected that the altitude where the bulge shifts hemispheres will increase. 

5. The high latitude winter helium bulge is expected to be caused by a 
seasonal variation of helium concentration in the lower thermosphere which 
extends into the exosphere . This seasonal variation may be related to a slight 
shift in the effective altitude at which diffusive separation of helium begins. 

6. The atomic oxygen bulge is expected to closely correspond to peak exo- 
spheric temperature . 
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Figure 2.- Relationship between observed values of parameter B, which gives latitudinal- 
seasonal variation of diurnal bulge position, and relative concentrations of helium to 
atomic oxygen. 
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